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Spin-dependent quark densities, matrix elements of specific density operators in 
proton states of definite spin-polarization, indicate that the nucleon may harbor 
an infinite variety of non-spherical shapes. We show that these matrix elements 
are closely related to specific transverse momentum dependent parton distributions 
accessible in the angular dependence of the semi-inclusive processes ep — > eirX and 
the Drell-Yan reaction pp — ► UX. New measurements or analyses would allow the 
direct exhibition of the non-spherical nature of the proton. 



Since the discovery that that the spins of quarks and anti-quarks account for only about 
33% of the nucleon spin [lj], 0], many experiments have sought the origins of the remain- 
der, which must be accounted for by effects of quark and gluon angular momentum. The 
importance of orbital angular momentum is also demonstrated in exclusive reactions. Mea- 
surements ji, H find that the ratio of the proton's electric and magnetic form factor Ge/Gm, 
falls with increasing momentum transfer Q 2 for 1< Q 2 < 6 GeV 2 . This striking behavior 
indicates that the sum of the orbital angular momentum of the quarks in the proton is 

non-vanishing a a a a. 

One expects that the presence of significant orbital angular momentum would lead to 
a non-spherical shape, if such can be defined by an appropriate operator. We showed [§], 
using the proton model of Ref. that the rest-frame ground-state matrix elements of 
spin-dependent density operators reveal a host of non-spherical shapes. The use of the 
spin- dependent density operator allows the detailed connection between orbital, spin and 
total angular momentum to be revealed in quantum systems. In the model of [10|| the 
orbital angular momentum originates from the relativistic nature of the quarks manifest by 
lower components of Dirac spinors of the wave function, but there are many other potential 
sources. 

It is natural to ask if the non-spherical nucleonic shapes can be measured. While ma- 
trix elements of the non-relativistic spin-density operator have been measured in condensed 
matter systems llj to reveal highly non-spherical densities and the related the orbital angu- 
lar momentum content of electron orbitals, finding the corresponding determination of the 
nucleon properties has remained a challenge. Our purpose here is to show that matrix ele- 
ments of the spin-dependent density are closely related to specific unintegrated (transverse 
momentum dependent) parton densities that could be obtained by measuring the angular 
dependence of the ep — > enX reaction and of the Drell-Yan production cross section in pp 
collisions. 

We begin by explaining how the shapes of a nucleon are exhibited by studying the rest- 
frame ground-state matrix elements of spin-dependent density operators [9J. The usual 
density operator in non-relativistic quantum mechanics is given by 

p(r)=5>(r-r 4 ), (1) 
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where is the position operator of the Vth particle. Matrix elements of this operator yield 
the density of a system. Suppose the particles also have spin 1/2. Then one can measure 
the probability that particle is at a given position r and has a spin in an arbitrary, fixed 
direction specified by a unit vector n. The spin projection operator is (1 + a • n)/2, so the 
spin- dependent density operator is 

p(r,n) = ^5(r-r J )i(l + CTi -n). (2) 

i / 

The spin-dependent density allows the presence of the orbital angular momentum to be 
revealed in the shape of the computed density. To understand this, it is worthwhile to 
consider a simple example of a single charged particle moving in a fixed rotationally invariant 
potential in an energy eigenstate of quantum numbers: I = l,j = 1/2, polarized in the 
direction s and radial wave function R(r). We find 

R 2 (r) 

p(r,n) = (W |p(r,n)| = — + 2cr • r n ■ r - cr ■ n|s). (3) 

Suppose n is either parallel or anti-parallel to the direction of the proton angular momentum 
defined by the vector s. The direction of the vector s defines an axis (the "z-axis"), and 
the direction of vectors can be represented in terms of this axis: s • r = cosd. With this 
notation p(r, n = s) = R 2 {r) cos 2 9, p(r, n = — s) = R 2 {r) sin 2 6 and the non-spherical shape 
is exhibited. The average of these two cases is a spherical shape, as is the average over the 
direction of s or the average over the direction of n. 

The densities of Eq. flTJ and Eq. d2J) can be extended to include other operators. Indeed, 



Ref. [ll( uses only the spin-dependent term appearing in Eq. (j2J), and this is weighted by the 
electronic charge. For systems of quarks, the densities could be weighted by the charge of the 
quarks, or be concerned with a specific flavor. One could also weight the spin-dependence 
by other operators. In particular, consider 

PR EL (r,n) = ^5(r- ri )^(l+7°^-n), (4) 

where the relativistic aspects are emphasized by the appearance of the Dirac operator 
7°, which becomes unity in the non-relativistic limit. We denote the density operators 
of Eq. (pQ),Eq. ([2]) and Eq. (J4]), and any number of obvious extensions, simply as densities. 

These densities are defined in terms of position, but to use QCD it is necessary to define 
operators that give the probability for a particle to have a given momentum, K, and a 
given direction of spin, n. The field-theoretic version of the spin-dependent charge density 
operator, Eq. (j2J), is a generalization of the operator defined in Ref. [§]: 

p(K,n) = / ^ e -^(0)( 7 ° + 7-n 7 5)£((U; path)^^)^, (5) 

where i/j is a quark field operator and flavor indices are omitted. The quark field operators 
are evaluated at equal time and accompanied by a path-ordered exponential link operator 
£((0, £; path) needed for color-gauge invariance. This introduces a path-dependence, which 
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must be specified correctly to obtain a parton interpretation. We will use the choice of 



Ref. [12|; see below. The first quantized version of Eq. §5§ (neglecting the gluonic aspects) 
is of the form of Eq. ([2]) except that the factor S(r — rj) is replaced by <5(K — Kj), where Kj 
is the momentum of the i'th quark. It is worthwhile to define another density corresponding 
to /5rel of Eq. (J4]) . This is given by 

p REL (K,n) = J ^ e -* K ^(0)7 (l + 7-n 7 5)A(U; path)^)^ . (6) 

The matrix element of a density operator in a nucleon state \P, S) of definite total angular 
momentum defined by four-vector and momentum P is 

p(K,n,S) = (P,S\p G (K,n)\P,S), p REL (K,n,S) = (P, S\p REL (K,n)\P, S). (7) 

The most general shape of the proton in its rest frame, obtained if parity and rotational 
invariance are upheld is then [9| 



p(K, n, S) = A(K 2 ) + B(K 2 )n ■ S + C(K 2 ) (n • K S ■ K - ■ S 
p REL (K, n, S) = A REL (K 2 ) + 5 REL (K 2 )n ■ S + C REL (K 2 ) (niS-K-^S 

with the last terms generating the non-spherical shape. Any wave function that yields a non- 
zero value of the coefficient C(K 2 ) or C RE l(K 2 ) represents a system of a non-spherical shape. 
If the relativistic constituent quark model of [10| is used, the principle difference between 
p(K, n, S) and pc(K,n, S) would be that C RE l = —C. This indicates that either C RE l or 
C can be used to infer information about the possible shapes of the nucleon. Measuring 
either C(K 2 ) or C RE l(K 2 ), would require controlling the three different vectors n, S and K 
or their equivalent. 

The densities of Eq. flHJ) are difficult to measure because the system must be probed 
without momentum transfer and the initial and final states are the same. This configuration 
also appears in parton distributions, both ordinary and transverse-momentum-dependent 
TMD. Parton density operators depend on pairs of quark-field operators defined at a fixed 
light cone time £ + = £ 3 +£° = while the density operators of Eq. (jSJ) and Eq. ([6]) are defined 
as an equal-time, £° = 0, correlation functions and cannot be regarded as parton density 
operators. However, we find a relation between the two sets of operators by integrating 
Eqs. f l5ll6l) over all values of K z . This sets £ 3 = 0, so the quark field operators of Eq. (j3j) or 
Eq. (jSJ) are now evaluated at £° = and £ 3 = so = 0: 

p r (K r ,n) = J^dK z p(K,n) = J |^ e - iK -^(0)( 7 ° + 7 ■ n 75 )£(0, £; n~)^ T ) ^ ±=q 
p RELT (K T ,n)= J^dK z p REL (K,n) = J |^ e - K -^ ^(0)t°(1 + 7 ■ n T5 )£(0, & n~)^ T )\^ ± j9) 
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where the specific path n~ is that of Appendix B of [121 ] . To obtain the relevant transverse 
densities we take the matrix element of par i n a nucleon state polarized in the transverse 
direction St, with 



p T (K T , n, S T ) = (P, S t |pt(K t , n) \P, S T ) 



(10) 



A T (K T ) + B T (K T )n ■ S T + C T {K T ) 



n ■ K T S T • K T - \K T n ■ S T 

M2 



PreL t(Kt, n, S T ) = A RELT (K T ) + B RELT (K T )n ■ S T + Crelt^I) ~ ' ' ^ ^^-f 



M 2 

where M is the nucleon mass, and we take the unit vector n to be in the transverse direction. 

Information regarding the shape of the nucleon resides in the functions CV, Crelt- We 
now are able to connect our newly-defined transverse densities with TMDparton distribu- 
tions. The latter are related to Dirac projections of correlation functions 12 ]: 



$ [r] (x,Kn 



2 (2tt) ; 



(12) 



£+=0 



The projections $^ depend on the fractional momentum x = K + / P + , Ky and on the hadron 
momentum P (in essence only P + and M where we work in a frame in which P + ^> M.). 
Depending on the Lorentz structure of the Dirac matrix T the projections <J>' r ^ are ordered 
according to powers of Mj P + multiplied with a function depending only on x and K T . Each 
factor M/P + leads to a suppression by a power in cross sections [13[ so that one may refer 
to the projections as having a 'twist' t related to the power [M / P + ) l ~ 2 that appears. Then 
moments (in x) of the KT-integrated functions involve local operators of twist t 14 . 

We use certain transverse momentum dependent parton distribution functions of 12) . 
which for transversely polarized nucleons are given by 



<f>^(x, K T ) = f 1 (x,K 



^ i+ ^(x, K T ) = S T h t (x, K T ) + ( KtKt ^ t ^) S ° t h ± AXy K 2 t)) 



(13) 
(14) 



$[^1(x,Kt) = ^9t(x,K t ) + g- l M2 } 9^K T ). (15) 

Terms of higher order in {Mj P + ) are neglected in the extraction of the functions gx, g? from 
high energy data. Similarly, at high energies, we may replace 7 + by v^T - Note that the 
quantity g? is closely related to the quantity CV, while the quantity h^ T is closely related 
to Crelt because ia t+r y 5 = 7+7*7 5 — > v / 27°7*7 5 - Extracting g^ would require a higher 
twist analysis, while hjj, appears at leading order in the cross sections for semi-inclusive 



leptoproduction experiments [15J. Thus the relativistic spin-dependent density Eq. ( ITTi) is 
easier to measure than the quantity of Eq. ffTOl) . 

We integrate the above parton distribution functions over all x so that the field operators 
are evaluated at ^ = as in our spin- dependent densities. A tilde is placed over a given 
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quantity to define the x-integrated result, e.g. $[ r l(Kj-) = / da;$' r '(i, Ky), fi(K T ) = 
J dxf\(x, Kip), etc. Then 



^ [rl(KT) = / e ~ lKT ^ T (P>s\mrc{o,bnJ)f/,(t T )\p,s) 



. (16) 

5+=o,?-=o 



The various $t r ' are expressed, in the infinite momentum frame, in terms of transverse 
momentum dependent parton distribution functions Eq. (]14|) and Eq. (fT5|) . To relate these 
to our functions At, Bp, Cp we evaluate these equations in the rest frame (P + — > Mj\pl, 
7 + — > "s/27 ) and choose the operators T to correspond to those appearing in the spin- 
dependent densities. We find that 



V2p(K T , n, S r ) = f x {K T ) + ~g T (K T )n ■ S T + ^ — - ^t(^I), (17) 



M 2 

1 ZV2, 



v^p RELr (K T , n, Sr) = h{K T ) + h x {K T )n ■ S T + • K r S r ■ Kr - ^■S T ) -, (f ^ 



Finding a non-zero value of either qt or /i^ T would demonstrate that the proton is not 
spherical. Both the spectator model 16| and the quark model T(| yield the result that 



g T = —h^ T . This relation does not appear to be a general result. In particular, one may 



use the most general parameterization of the correlation functions [121 . |17| in terms of scalar 
functions Ai(K, P) to show that g^ is proportional to the term As and hyj. is proportional to 
the term P + / M An. This means that, in the nucleon rest frame, the ratio g^ /hyp is simply a 
function of (x, Ky). Furthermore, the quantity h^ T is known to characterize the dependence 
of the transverse polarization of quarks in a transversely polarized nucleon on the direction 
of k t m, in . Thus the two quantities g^ and hj-j, each characterize the non-spherical 
nature of the nucleon, and the density prelt can be thought of as "the" spin-dependent 
density" . 

Next we focus on experimental means to access hyr- The presence of the term hyp causes 
distinctive signatures in semi-inclusive leptoproduction experiments [l5| in which a hadron h 
is produced. If the target is polarized in a direction transverse to the lepton scattering plane, 
the cross section acquires a term proportional to cos(3(f) l h ) where <j) l h is the angle between 
the hadron production plane (defined by the momenta of the incoming virtual photon and 
the outgoing hadron) and the lepton scattering plane. A similar effect occurs in electroweak 
semi-inclusive deep inelastic le ptop roduction and this could be accessed at high energies 
such as those found at HERA [20|]. Another signature occurs in the angular distribution 
of the leptoproduction of p mesons (2l|], obtained using an unpolarized lepton beam and 
a transversely polarized target. Similarly the term hyp makes its presence felt in studying 
the production of two-pions inside the same current jet (22|. In each of these cases, the 
momentum of the virtual photon and its vector nature provide the analogue of two of the 
three vectors n and Sr needed to define the spin-dependent density. The hadronic transverse 
momentum provides the third, K r . 



FIG. 1: (Color online) Transverse shapes of the nucleoli: y/2pr('K.T,' a )/ fi(K T ). The horizontal 
axis is the the direction of St and n = St, <fi n = 0. The shapes vary from circular to highly 
deformed as Kt is increased from to 2.0 GeV in steps of 0.25 GeV. 



FIG. 2: (Color online) Transverse shapes of the nucleon, as in Fig. 1 except that <j) n = ir. 

Another interesting possibility occurs in the Drell-Yan reaction pp{\) — > MX using one 
transversely polarized proton [23[. In case the term hj- T causes a distinctive oscillatory 
dependence on the angle 30 — (j>s 1 , where <fi is the angle between the momentum of the 
outgoing lepton and the reaction plane in the lepton center of mass frame, and <ps 1 denotes 
the direction of polarization with respect to the reaction plane. The term h^ T is multiplied 
by the anti-quark Boer-Mulders function. 

To illustrate the shapes that could be obtainable using this method we use the spectator 
model of [l6| to evaluate the shapes of the proton. We rewrite Eq. ffl8|) as 



v / 2pRELT(K r ,n) Mgjj 1K T Mt(K t ) hQ , 

~ 7. = 1 + — — COS0 n H — — COSIZ0 — On)— — , (19) 

h{K T ) fi(K T ) 2M^ yv ' fx{K T ) K ' 

where is the angle between and St and <p n is the angle between n and St- The 
transverse shapes of the nucleon, as defined by the right hand side of Eq. (1191) are shown in 
Fig. 1, taking <j) n = 0. Deformation is seen for values of Kt as small as 0.25 GeV, and this 
increases as Kt increases. Choosing <fr n = ir emphasizes the non-spherical nature because 
the first two terms of Eq. ( Tl9l) tend to cancel. The possible shapes implied by Eq. (Tl9l 
can be thought of as transverse projections of the shapes displayed in Refs. [9jj. One 
complication as that it would be very difficult to measure the necessary TMD's at all values 
of x to construct the integrals appearing here. However, the model [l6|] indicates that the 
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functions fi, hi and h± T have very similar x dependence, so that measurements at values of 
x for which these functions peak should be sufficient. 

We have shown that that the non-spherical nature of the nucleon shape is closely re- 
lated to the non-vanishing of the measurable TMD hj~ T . While determining this function 
experimentally represents a challenge, the ultimate determination of a non-zero values would 
clearly demonstrate that the shape of the proton not round. 
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